ASYMPTOTIC STATISTICAL CHARACTERIZATIONS OF p-HARMONIC 
FUNCTIONS OF TWO VARIABLES 
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Abstract. Generalizing the well-known mean-value property of harmonic functions, we prove that 
a p-harmonic function of two variables satisfies, in a viscosity sense, two asymptotic formulas 
involving its local statistics. Moreover, we show that these asymptotic formulas characterize p- 
harmonic functions when 1 < p < oo. An example demonstrates that, in general, these formulas do 
not hold in a non-asymptotic sense. 



1. Introduction 



< 

A fundamental and fascinating fact about harmonic functions is their characterization by the 
mean value property 01 : the continuous function u is harmonic in the domain O. c R N if and only 
if 



u(x) = -r u(s)ds = t u(y)dy for each xeQ, (1) 

JdB r (x) JbAx) 



where B r (x) € Q is a ball with center x and radius r > 0, dB r (x) is its boundary, and JF / denotes 
the average of / over the set E. Ostensibly, identity ([I]) says nothing about derivatives and could 
be studied entirely within the category of continuous functions. It is the prototypical statistical 
characterization of solutions of a PDE, and it is natural to wonder if this is peculiar to Laplace's 
equation. In other words, can one characterize solutions of other PDEs in a statistical way that 
avoids any explicit mention of derivatives? 

Recent work shows that such statistical characterizations exist, in a certain sense, for p-harmonic 
functions, i.e., solutions of the quasilinear PDE 



— A p u : = — div (\Du\ p - 2 Du) = , for 1 < p < oo. (2) 

More precisely, p-harmonic functions are usually defined to be weak solutions of ©; thanks to 
work by Juutinen et al. (HI, however, weak solutions of © are the same as viscosity solutions of 
©. Viscosity techniques are particularly relevant to the present work, as Manfredi et al. [10J used 
such methods to prove that the continuous function u is /^-harmonic in the domain Q c M, N if and 
only if the functional equation 

u(x) = %■ \ max u + min u \ + B t u(y) dy + o(s 2 ) as s — > (3) 

2 \B s (x) B s (x) j JBAx) 
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holds in the viscosity sense for all x e Q. The constants a and f3 are determined by the exponent p 
and the dimension N: 

p-2 , „ 2+N 

a := and yS := . 

p + A/ /? + A/ 

This characterization also holds for oo-harmonic functions, where the oo-Laplacian A^ has the 

formal definition 

^ 1 yi du du d 2 u 

\Du\ 2 dxi dxj dxidxj 

for smooth u. 

To establish their results, the authors of IfTOl combine several interesting facts. First, calculating 
formally yields 

A p u = \Du\ p ~ 2 (Aw + (p - 2)A 0O ) , (5) 
an identity that plays a central role in both (8]| and IfTOl . Using it, Juutinen et al. proved that u is a 
viscosity solution of © if and only if 

-Am - (p - 2)A co u = 

in the viscosity sense, about which more will be said below. Manfredi et al. then invoke the 
identities 



2 



u(x) - f u(y) dy = - — — Au(x) + o(s 2 ) (6) 

JB E (x) z vV + l ) 

and 

If ) 2 

u{x) - — \ max u(y) + mm u(y) > = -—A^uix) + o(e ) , (7) 

2 [yeBAx) ' yeB s (x) ) 2 

valid for smooth functions as e —* 0, to obtain their asymptotic characterization ©. Here and in 
what follows, a function is called smooth if it is C 2 . 

The decomposition © can be written in various ways, a fact that we exploit to obtain new 
statistical characterizations of p-harmonic functions of two variables. Specifically, if we define 
the 1-Laplacian A t on smooth functions by 

/ Du \ 

A 1 »:=|DMd W ( lsr ). (8) 

then the formal relationship 

Ai = A - Aco 

holds and immediately yields two alternatives to © : 

A p u = \Du\ p - 2 {{p-l)Au + {2- p)A lU ) , (9) 

and 

A p u = \Du\ p ~ 2 ( AiU + (p — l)AooW) . (10) 
Using these identities and the Taylor approximation 

s 2 

u(x) - median {u(s)} = --Aiu(x) + o(s 2 ) , (1 1) 

s&WAx) 2 

valid for smooth functions u of two variables as e — » 0, we prove the following: 

Theorem 1. Suppose that 1 < p < oo and OcR 2 is open, and let u be a continuous function on 
Q. The following are equivalent: 
(1) u is p-harmonic in Q. 
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(2) At each x e Q, the equation 

u(x) = (— - 1 ) median} } + |2 - — ) -f u(s)ds + o(s 2 ) as e — > (12) 

holds in the viscosity sense. 

(3) At each x e Q, f/ze equation 

u{x) = — median { u(s) } + | ^- - 1 1 max { u(y) } + min { w(y) } ) + o(e 2 ) as e — > (13) 

^ ,e^ £ (x) \ 2p )\yeBM ycBM J 

holds in the viscosity sense. 

The median operator occurring here is defined as expected: if u is continuous on Q., x e Q, and 
B £ (x) c n, 

if and only if 

| {s e dB E {x) : u(s) >m}\ = \{se dB E {x) : u{s) < m } | , 

where |£1 is the 1 -dimensional Hausdorff measure of the set E. We remark that if u is smooth and 
\Du(x)\ ± 0, then (fl"2l) and (fl"3l) hold in the usual non-viscosity sense if and only if A p u(x) = 0. 
This follows from Lemmas 1 and 2 below. 

Considering (QQ), it is natural to ask if the formulas (fl"2l) and (fl"3l) hold in a non-asymptotic sense. 
More precisely, if u is p-harmonic in Q., do the equations 

u{x) = | - - 1 ) median { w(^) } + |2 - -J -f u{s)ds (14) 

PI JdB E (x) 



m = median { u(s) } 

sedB £ (x) 



p ) sedBAx) 



u(x) = — median { u{s) } + \ ^ \\ max { uiy) } + min { u(y) } I (15) 

p sedB c (x) \ 2p j \yeB s (x) yeB e (x) 

necessarily hold at all x e Q, for all s > sufficiently small? The answer to this question is no, and 
in Section 12.31 we provide an example demonstrating that these equations do not hold in general 
even for smooth p-harmonic functions. 

On the way to proving Theorem \T\ in Section 12.21 we provide a simple analytic proof of iden- 
tity (fTTT) . We should point out, however, that the relationship between median values and the 
1-Laplacian has appeared before, either explicitly or implicitly. In [12], for example, Oberman 
uses a discrete median scheme of forward Euler type to approximate solutions of the parabolic 
mean curvature equation, 

du 

- — Aim = for t>0, u(-,0) = u , (16) 
ot 

in two space dimensions. Unlike many other proposed algorithms for this equation, Oberman's 
median scheme is provably convergent, an easy consequence of the main theorem in 0J. 

Kohn and Serfaty [9] discuss a different convergent approximation scheme for the initial-value 
problem ([TBI that can be described geometrically as follows. Let T(0) be a simple closed curve in 
the plane, let Y(t) be the curve obtained from T(0) by letting it evolve by mean curvature for time 

2 

t, and fix a small s > 0. The curve F(t + y) is approximately the locus of all centers of circles 

2 

of radius s with antipodal points on T(t); one can approximate Y(t + — ) by tracking the center of 
a segment of length 2s as its endpoints traverse the curve T(t). This is the basic idea behind our 
proof of (fTTT) . even though Kohn and Serfaty never mention medians in Q- Related papers that 
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use similar ideas without explicitly connecting the 1-Laplacian and median values include, but are 
certainly not limited to, @ and IfTTTl . 

The present work is actually closely related to the work of Jackson and it is our pleasure to 
briefly discuss this connection. Over the past thirty or so years, viscosity solutions have become 
a standard tool in the study of nonlinear PDEs. However the contemporary viscosity approach is 
similar in some ways to the earlier abstract Perron method of Jackson and Jackson and Beckenbach 
as in J6l and 0. In fact, for a class of second-order elliptic PDEs, viscosity subsolutions 
and the subfunctions of Beckenbach and Jackson are equivalent (see [5]). Furthermore, Jackson 
applied this abstract Perron method to obtain existence and uniqueness results for the minimal 
surface equation in two independent variables Q; this work is closely related to ongoing work 
on 1-harmonic functions IfTBl , as the level sets of 1-harmonic functions are minimal surfaces (cf. 
d). 

2. New results 

2.1. Definitions. Before proving Theorem [H we review the necessary definitions and related re- 
sults. 

Definition 1. Suppose that 1 < p < oo, and let Q be a domain in K 2 . 

(1) The lower semicontinuous function u is p-superharmonic in CI in the viscosity sense if and 
only if the equivalent inequalities 

(1 - p)A(p + (p - 2)A 1 cp > and - A Y (p + (1 - p)A DO cp > (17) 

hold at x e CI for any smooth function </? such that \D(p(x)\ ^ and u - ip has a strict 
minimum at x. 

(2) The upper semicontinuous function u is p-subharmonic in CI in the viscosity sense if and 
only if the equivalent inequalities 

(1 - p)A<p + (p - 2)A l (p < and - A\<p + (1 - p)A co (p < (18) 

hold at x e CI for any smooth function <p such that \Dtp{x)\ + and u - ip has a strict 
maximum at x. 

(3) u is p-harmonic in CI if it is both p-superharmonic and p-subharmonic in CI. 

The legitimacy of this definition follows from JH and the formal identities ©, ® and (flOl) 
above, as checking p-harmonicity in the viscosity sense reduces to evaluating -A p (p for smooth 
functions (p away from critical points. We refer to [HI and IfTOl for more details. 

Definition 2. Let 1 < p < oo, let CI be a domain in K 2 , and consider the equation 

u(x) = \- - 1 1 median {«(■?)} + ( 2 - -) -f u(s)ds + o(s 2 ) as s -> 0. (19) 

\p ) S €dB B (x) \ pJJdBAx) 

(1) u is a supersolution o f/[T9\) in the viscosity sense if and only if the inequality 

(fix) >{-- 1 ) median { (p(s) } + ( 2 - - 1 £ <p(s) ds + o(e 2 ) as s (20) 
\p ) sedBAx) \ p)J dBs ( X ) 

holds at x 6 CI for any smooth function <p such that \D(p(x)\ ± and u - cp has a strict 
minimum at x. 
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(2) u is a subsolution of([[9!i) in the viscosity sense if and only if the inequality 

(fix) <{-- 1 1 median { (p(s) } + (2 - - ] -f w(s) ds + o(s 2 ) as £ (21) 

\p ) sedBAx) \ P) J d B E (x) 

holds at x 6 Q, for any smooth function (p such that \Dip(x)\ ^ and u - tp has a strict 
maximum at x. 

(3) u is a solution of 4791) z>z ?/ze viscosity sense if and only if it is both a subsolution and a 
supersolution. 

2.2. Proof of Theorem [IJ We begin with asymptotic formulas valid for smooth functions that 
will be used to establish our main result. The following lemma can be established using Taylor 
expansion; we omit the routine proof. 

Lemma 1. Let Qbe a domain in R 2 , let jefi, and let <p be a smooth function on Q. Then 



s 2 



(p(x) - -r ip(s) ds = -—A(p(x) + o(e z ) as s 0. (22) 

JdBAx) 



Lemma 2. Let Vibe a domain in R 2 , let x = (xi,X2) 6 Q, and let <p be a smooth function on Q. with 
\Dcp(x)\ * 0. Then 

s 2 , 

(fix) - median {w{s)\ = — — Aiipix) + o(s ) as e — > 0. (23) 

sedBAx) 2 

Proof. The Implicit Function Theorem guarantees that, for e > sufficiently small, the level sets 
of (p form a one-parameter family of smooth, non-intersecting curves that foliate the closed ball 
B £ (x). Consequently, the median of <p over dB £ (x), 

M E : = median {w(s)} , 

sedB s (x) 

is the value corresponding to the level set that intersects dB E {x) in antipodal points; for each s > 0, 
there is a unique angle 6 e e [0, 2n) such that 

M E = (p(x\ + s cos 6 E , x 2 + s sin 8 £ ) = y{x\ - s cos 6 E , x 2 - s sin 6 E ) . (24) 

Let v £ denote the unit vector (cos 8 £ , sin 6 £ ), and define 

Dip\x) :=(-^ 2 (x),^(x)). 

The derivatives of (p below are evaluated at x, which we omit for simplicity. Taylor expanding 
about x yields 



e 2 

M £ = <p(x + s\ £ ) = cp(x) + sD(p ■ \ £ + —y T £ D 2 <p\ £ + o(s 2 ) (25) 



and 



s 2 

M £ = (p(x - s\ £ ) = (p(x) - sD<p ■ \ £ + —\lD 2 (p\ £ + o(s 2 ) . (26) 



Since these expressions both equal M £ 
We therefore have 



sD(p ■ \ £ = o{s 2 ) . 



|Zty| 
where 

sD(p ■ w £ = o(s 2 ) , 



V +w £ , (27) 
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and we see (among other things) that the sequence {v e } of unit vectors converges: 

v £ "77TT as £ J- • 
\Dip\ 

Using the decomposition (12"7T ) in the right-hand side of either (1231 ) or (|26l ) yields (cf. @) 



£ 2 (D<p- L y 2 2 e 



2 



p(x) - M e = ~^-f - D z (f -Z- + o(e z ) = -—Akp + o{s 2 ) , (28) 
proving the lemma. 

□ 

With these lemmas, Theorem \T\ is easily established using the same approach as in IfTOll : apply 
the asymptotic formulas for smooth functions to the viscosity formulation. 

Proof. Suppose that u is continuous in Q. and that (p is a smooth function for which \D(p(x)\ ± 
and u - ip has a strict minimum at x e Q.. Using Lemmas 1 and 2 and observing that (2/p - 1) + 
(2 - 2/p) = 1, it follows that the first inequality in (TT71) holds if and only if ([2~0l holds. Thus u is 
p-superharmonic in the viscosity sense if and only if it is a viscosity supersolution of (fT21 . The 
analogous argument establishes the equivalence of p-subharmonicity and being a subsolution of 
CU). 

The equivalence of the first and third statements of the theorem is proved similarly, using identity 
CO) instead of Lemma 1 . 

□ 

2.3. Necessity of Asymptotic Nature of Theorem Q] In this section, we present an example to 
show that (fT4l) and (TT5T) do not hold for p-harmonic functions in general. In fact, these equations do 
not even necessarily hold for all s > sufficiently small, so that the asymptotic results appearing 
in Theorem [T]are, in general, the best available. 

For any 1 < p < 2, the function u p (x) = \x\ < - p ~ 2)l<p ~ 1) is smooth and p-harmonic in R 2 \ {0}, and is 
known as the fundamental solution of the p-Laplacian (see for example [El). Let x = (x\, 0) where 
X\ > and let < e < X\ . Because u p is radial and radially decreasing, it is not hard to see that 

median u„ = (x\ + £ 2 )(p- 2 ^p~ 1 ) . (29) 

dB s (x) 



The mean of u p on dB E (x) is 

^-j (x] + 2XJSCOS 6 + e 2 )^- 2 )/ 2 ^- 1 ) de. (30) 
Using (|29l ) and (1301) . (fl4l) at * with u = u p becomes 

\ Xl \WW = 1^- lj^+e 2 )^^ 1 ^^ - i- ^(x 2 +2x ie cos0+ e 2 ) ( ''" 2)/2( ^ 1) (31) 



If (T3Tl) holds for all £ sufficiently small we can differentiate it with respect to e to obtain 

(2-p)(x 2 l +s 2 ) ip - 2)/2ip - 1) - l £=^--^ [ {x 2 l +2x l sco%e + s 2 ) {p ~ 2)l2(p ~ l) -\x l co^e + s)d6. (32) 

27T Jo 

Now let x\ = 1 and p = 3/2. The last equation is then 

-1 f 2?r 

(1/2)(1 +e 2 )" 3/2 e= — (1 + 2s cos 6 + s 2 T 3 12 (cos 6 + £)de, (33) 
2?r Jo 
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which holds if and only if 

1 r 2n (l + 2ecosfl + e 2 \~ 3/2 1 C 2n I 2s cos OY 3 ' 2 

-£=- - (cos6 + s)d6= - \ 1 + - (cos6 + s)d6. (34) 

Jo \ 1 + e 2 / Jo \ 1 + e 2 / 

Using the binomial formula: 

, 2ecos#\" 3/2 , 3/2ecos#\ 15 /2ecos0\ 2 35 (2ecos6\ 3 

1 + ^^ =1-X T-T ^-^T -TT T7T (35) 



1 + e 2 / 2 \ 1 + e 2 / 8 \ 1 + e 2 / 16 \ 1 + s 

plus higher order terms. Therefore the integrand in (1341) is equal to 

3e 2 15e 2 cos 3 # 35e 3 cos 4 # 3e 2 cos# 15£ 3 cos 2 # 

cos ,___ cos + (36 ) 

plus terms of order 4 and higher. Using (1361 ) in the integral in (1341) , noting that odd powers of cos 6 
integrate to zero and recalling that JC cos 2 6 d6 = n and L cos 4 6 d9 = (3/4)n, we obtain 

1 r 2n (l+2scos9 + e 2 \~ 3/2 , 

(cos0 + e)J0« -e-(21/8)e 3 , (37) 

tJo \ 1+e 2 / 

which is strictly less than —e if e is sufficiently small so that (1341) does not hold. As a result, (T311) 
cannot hold for all s sufficiently small. 

The same example can be used to show that (fl"5l) also fails in general, even if e is small. Again 



let p = 3/2 and x = (1,0), and let < s < 1. The maximum value of w p on 5 £ (jc) is 1/(1 — e) and 
the minimum on the same ball is 1/(1 + s). Using (|29l) , in this case (fT5l) becomes 



3 V / 6\l-e 1+e/ 



(38) 



which one can easily see does not hold, even if s > is restricted to being smaller than some e . 



3. Concluding remarks 

The asymptotic characterizations of p-harmonic functions in [ T0j| are valid in N dimensions. It 
would be interesting to extend the results presented here to higher dimensions. The only part of 
the proof of Theorem 1 that requires two dimensions is Lemma 2. If an A/-dimensional version of 
Lemma 2, perhaps involving the median on an (N - l)-dimensional sphere, were established, new 
asymptotic statistical characterizations of p-harmonic functions would follow. 

We presented an example showing that, in general, only asymptotic characterizations of this 
type are possible. However, this is not the case for p = 2. A natural question is: do the equations 
(fT4l) and (fl5l) hold either globally or locally for any other values of pi Concrete examples in [fT3l 
show that the limiting cases of (fT4"l) and (fl5l) can hold when p = 1 , but more work on this question 
needs to be done. 

Finally, we did not consider the extreme cases p = 1 and p = oo, although we remark that if p is 
formally allowed to be oo in (fTSl the resulting characterization is the same as that in [[101 . 
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